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Abstract. We aim at proving existence of weak solutions to the stationary compressible Navier- 
Stokes system coupled with the Allen-Cahn equation. The model is studied in a bounded three 
dimensional domain with slip boundary conditions for the momentum equations and the Neumann 
condition for the Allen-Cahn model. The main result establishes existence of weak solutions with 
bounded densities. The construction is possible assuming sufficiently large value of the heat capacity 
ratio J ip ^ q'^)- As a corollary we obtain weak solutions for a less restrictive case losing point-wise 
boundedness of the density. 


1 Introduction 


Phase transition phenomena are important subject of applied mathematics. Complexity of physical 
nature of these processes makes us to choose different models in order to obtain the best description 
in a concrete studied case. Here we want to concentrate our attention on phenomena with fuzzy 
phase interfaces. We think about the process of melting or freezing, at the level of almost constant 
critical temperature. We want to control densities of phase constituents, and here we arrive at the 
Allen-Cahn equation: 


df, . . 

Q^iQ.c) - Ac = Qp. 


( 1 ) 


The equation above gives us control of one of the phase constituents in terms of the chemical 
potential p. In the chosen model this equation is coupled with the compressible Navier-Stokes 
system, which represents a well established and frequently studied model describing the flows of 
viscous compressible single constituted fluids. 

The mathematical analysis of coupled systems: the compressible Navier-Stokes type and the 
phase separation is in its infancy [HEIITKIIIIITIIIB], although the mathematical theory of each of 
them separately is quite developed (see e.g. [9], mm)- In this article, we study existence of steady 
weak solutions. The thermodynamically consistent derivation of the model under consideration, 
which is a variant of a model proposed by Blesgen [2], was presented by Heida, Malek and Rajagopal 
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in m- It is represented by the following system of partial differential equations for three unknowns: 
density of the fluid g, velocity field v, and concentration of one selected constituent c 

div(pv) = 0, 

div(pv (g) v) = divT + pF, 

div(f)cv) = —fj., 

df, s . 

QH = g-^{Q,c) - Ac, 

where the stress tensor T is given by 

T(Vv, Vc, g, c) = §(Vv) - (Vc (g) Vc - - p{q, c)I, 

with the thermodynamical pressure p{g, <;) = l!% and the viscous stress S satisfying the Stokes 
law for Newtonian fluid 


( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 


S(Vv) = 1^ (Vv + V^v - div vl) + rj div vl, 

with the shear viscosity coefficient > 0 , and the bulk viscosity coefficient rj > 0 . The free energy 
is assumed in the form with the so-called logarithmic potential 

f{Q,c) = + (aic-F 02(1 - c)) log p-F clog c-F (1 - c) log(l - c) + b{c) 

with some 7 > lH, 01,02 > 0 and b a smooth bounded function with |6'(c)| < C. Moreover we 
assume without loss of generality that oi >02, and we denote for the sake of simplicity o = oi — 02, 
d = oi, and L(c) = clogc -|- (1 — c) log(l — c). The logarithmic terms related to the entropy of the 
system assure that c E [ 0 , 1] almost everywhere, since 

— (p,c) =logc-log(l-c)-h(oi - a2)\ogg + b’{c) 

=L'{c) -I- ologp -I- 6'(c), 
p{g, c) =g'^ + g[ac + d). 

The fluid is contained in a smooth bounded domain fl, we supply the equations in the domain 
with boundary conditions 


V • n = 0 at do., (6) 

n • T(c, Vv) • -|- fev • Tn = 0 at dVl, ( 7 ) 

Vc • n = 0 at (8) 

where parameter k > 0 represents the friction on the boundarjH, r^, n = 1,2 are two linearly 
independent tangent vectors to dkl, and n denotes the normal vector. 

^We could assume also more general pressure law function similarly to [ 10 ], but it would lead only to additional 
unnecessary technicalities, so we omit it. Similarly, we can consider more general singular functions L{c). 

^The values of 7 which ensure the existence result will be specified later. 

®For slip boundary condition corresponding to the case fc = 0, we need to assume that the domain ft is not axially 
symmetric. 
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The solutions are parametrized by means of the condition 

J gdx = M. (9) 

n 

The fluid is driven by an external force F E 

The main result of the paper is the following. 

Theorem 1. Let 7 > 6 , M > 0 and F E L°°(r2,R^). Then there exists at least one weak solution 
to the system such that c E [0,1] in n 

g E V E and c E for p < 00 . (10) 

In addition, if we assume only 7 > 3, then there exists at least one weak solution such that 

geL^^-^{n), vElT^’^(fl), VcEL^(fl). ( 11 ) 

The meaning of the solutions is explained in Definitions^ below. 

The existence of non-stationary solutions to the above mentioned system with no-slip boundary 
condition for the velocity was shown by Feireisl et al. in m- Ding et al. studied in [7] the global 
existence of weak, strong, and even classical solutions in ID with free energy approximated by 
a suitable bistable function, assuming no vacuum zones in the initial data. Further, Kotschote 
m put his attention to a more advanced model where the extra stress tensor is multiplied by 
density function. He showed, however only, the local-in-time existence of strong solutions provided 
positiveness of density, including the thermal effects as well. The existence of travelling waves for 
this model was shown by Freistiihler in m- To best our knowledge there are no results concerning 
existence of weak solutions to the stationary system ([2])-([5]). Concerning the steady solutions to 
the compressible Navier-Stokes system, we refer to the pioneering work of Lions [20], and further 
extensions [3lll2lll6ll27H29ll31] , we use a technique developed by Mucha, and Pokorny (see [22l[23| ), 
which was modified for the Navier-Stokes-Fourier system [24] . too. Methods from [22| allow to 
obtain solutions with point-wise bounded densities, which seems to be the best possible regularity 
for weak solutions. Here we find an explanation of titled decently regular solutions. 

We introduce the following definition of weak solutions to the system ([^- ([5]) . 

Definition 1. Let M > 0 be a given constant, 7 > 3, we say that quadruple g,v,p,c is a weak 
solution to the steady Navier-Stokes-Allen-Cahn system, if E L'*'(H), > 0 a.e. in H, v E iT^’^(H), 

p E L^(H), c E iT^’^(H), 0 < c < 1 a.e. on {g > 0}, with gL'{c) E (H) and v • n = 0 

satisfied on dLl in the sense of traces, and if the following holds true: 

1. Continuity equation is satisfied in the distributional sense, id est0 

div[g\r) = 0 , in P'(R^) 

2. For every ip E C°°(H,R^), (/? • n = 0 at dQ 

/ (~^’(v <81 v) : Vip + T(^, c, Vc, Vv) : V(p )dx -|- / k{'v ■ T){(p ■ T)dS' = gF ■ (p dx. 
n an u 


‘^Note that we do not need the notion of renormalized continuity equation, since our density is regular enough. 
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3. For every Lp G 


and 


/ 


- Vcipdx = / —fi(pdx 


/ 


J QpLip dx = J Q 


df{g,c) 
dc 


ip + Vc • Vipdx. 


( 12 ) 


such that f gdx = M. 

n 

The paper is constructed as follows. First we compute formal a priori estimate, it allows us 
to determine expected regularity of sought solutions. In Section [3] we deal with the approxima¬ 
tive system, we construct regular approximative solutions together with required estimates in the 
dependence of approximative parameters. Finally we analyze the limit, showing the strong con¬ 
vergence of approximative densities. The method is based on the fact that for the approximative 
densities we are able to find such m that 

lim > m}\ = 0. (13) 

e-S-O 

We skip the proof of the second part of Theorem [H we explain it in Remark [2j Through out the 
paper we try to follow the standard notation. 

2 A priori estimates 

Before the technical part of the proof, we will present here the a priori estimates on certain norms 
of the solution derived by purely heuristic approach. All generic constants, which may depend on 
the given data as well, will be denoted by C, its values can vary from line to line or even in the 
same formula. 

Lemma 1. Let Ll C be a bounded domain with boundary. Assume that all the above 

mentioned hypotheses are satisfied with 7 > 3 and that g,v,fi,c is a sufficiently smooth solution 
satisfying with the boundary conditions Then 

Il’'^llivl’2(r2) + ll£'lll,37-6(n) + ||Vc||^67-12 -|- \\h\\L'^(n) + \\qL' < C, (14) 

where C may depend on the data, but is independent of the solution. 

Remark 1. Note that > 2, for 7 > 3 and that from the bound of the last norm on the left 

hand side, we immediately deduce that c G [0,1] a.e. on the set > 0}. Moreover, if 7 > 4, then 
according to the Sobolev imbedding c is continuous function and we conclude from the maximum 
principle for harmonic functions that in fact c G [0,1] a.e. in Ll. 

Proof. First, multiplying the momentum equation by the velocity field v yields (with usage of 
Korn’s inequality and boundary condition) 

— J pdivvdx d- C J |Vvpdx < J gfiVc ■ v — g^Vc • v -|- ^>F • v dx, (15) 

SI SI SI 
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where we have used the equation ([5|) as well. Further, we conclude from the continuity equation 
that 

J pdivv + div(^3/v)dx = J g-^Vc-vdx, (16) 


o o 

and according to the constitutive equation for /r we get 

/ //^dx = / -QfiVc-vdx. 


(17) 


Thus, summing up (fT^ - (fT71) yields for p > 


/' 


7 


Vvf + p'^dx < C / • vdx < Cl 1 + 


'L5(0) 


<C 1 + 


I 3(p-l) \ 

\Lp{n) J ■ 


In order to bound the density by means of the Bogovskii estimates, we need Vc in L'^(ll) ,q>2. 
This can be deduced from the constitutive relation for p. We state this for purpose of future 
references more generally in the following lemma. 

Lemma 2. Suppose that f{Q,c) is as above, and p E L^{Q) and g E LP{Q) with q > 2, p > 3, 
q < satisfy the equation (0) with the boundary condition Then 


77 

'Sc 


(O) 




Lp^ (Q) 


Lip+Jq-pq (f2) 


Proof of Lemma\^ Let us test the corresponding equation by F(^L'{c)'j, for increasing function F 
with growth C-^(C) ~ with some /3 > 0, so (recall that L is convex) 

< I \Vcf F'{L\c))L"{c) + gL'{c)F{L'{c))dx = 
n 

= J g[p — a log g — b'{c))F(^L'{c))dx 

Q 

< Cy gFfi |L'(c)|^ (|/i| pTPr + gj+i |logp|) + g |L'(c)|^ |c|dx 


< C 


gP+^L'{c) 


L/3+1 


L. + ( / ^?|log^7+Mx)"^^) + p7^L'(c) 

L<J-(/3+l) \J 2 J 

n 


d+\ 

L/3+1 


which yields 




L/3+i{f2) 


< c 1 + 


LI + ( / £'|log£» 7 +Mx 

n 


1 

/3+1 


Now we fix the exponent (3 such that = p E (3, oo), hence 
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Thus, 


|oLTc)|| pi = 


^/ipq-g) g{pq-g-p)/ipg-<}) 


PQ 

Lp+p (n) 


< 

<C^I 

<CI 

<c 


gP/(P<}-<l)L'{c) 


g(p-l) 


Jpq-q-p)/{pg-q) 



L p (n) 

pq-p-q 

f ^ 

<j(p-i) 

Lp{n) 

+ A'' L9(o) 

pq-p-q 

/ p 

q(p-i) 

LP(Q) 

1 + 1,9(0) WsWlp^Q) 


pp(q-i) 
Lpi-p-1 (Q) 

(p-i-)Q 


I 9(P 




\Li(n} + 1 ) IkllLP(n) + 1 

For the other terms we have 

+ llellwcn))- 

- ll£’llLp(r2 ) II^IIl9(q) ■ 


Using the classical elliptic estimates on equation —Ac = together with embedding 


1 P'' 3pi7 

W ’p+i{Q,) L^p+^i-pi {Vl) we get 


||Vc|| 3vq < 

L3p+Jq-pq (57) 


df 

+ Q^J- 
oc 


qp 

L^(Q) 


+ l|c|l 7 ,s(Q) < C 


Li{n) 


+ 1) ll^llLp(r2) + ’ (1^) 


where for fl we can take {g > ^/o} which has positive measure and on which c € [ 0 , 1 ] a.e. according 
to the logarithmic terms in L'{c). This completes the proof of Lemma [2j □ 

Applying Lemma[2]with (7 = 2 yields ||Vc|| 6 p < C(l +. Now, we are ready 

L ‘s+P (O) 

to perform the Bogovskii estimates, id est to test the momentum equation bjo 

$ = -{(?“} 

where a > 0 will be specified later, and B ~ div“^ is the Bogovskii operator. The theory implies 
that ||V$||lp < C'||^>“||lp with = 0 (see e.g. [^[51127] 1 : we obtain 


Jp{Qj c)^>“dx < J (—(f?v ® v) : + S(Vv) : — ^>F • $)da 


n 


+ cl\Vc\^\Vt\d^ + J p{Q,c){Q'^}ndx. 


The terms on the left hand side of the inequality have sign and give the desired estimate of if 

the right hand side will be estimated, thus we set p = 7 + a. We will present only the most difficult 

®We use notation {< 7 }^ ~ '\n\ Iei9 
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and restrictive terms. 


/"iVcI^lV^ldx < ||Vcf 6p ||V$|| 3p 

J L^in) 

n 


<c(i + ||e|| 


(7p-6)/(3p-3) 

Lp{n) 


3p ) 


<c(l + ||£.|| 


(7p-6)/(3p-3) 

Lp{n) 


3otp—2p+6 , 


<c(l + ||(>ll 


(5p+3Q:p)/(3p—3) 
LP{n) 


(19) 


provided 2p-e — Pi equivalently 

0 < a < 2(7 - 3 ), ( 20 ) 

which yields the restriction 7 > 3 . The condition < 7 + a is satisfied even for all 7 > |, 

so we can put this term to the left hand side. 

Further, provided a < 27 — 3 


J (£»v( 8 ) v) : V$dx < ||v||^ 6 (n) II^372(^2) 

n 


<C{i + llell 


p 

3(P-1) 'i 

LP(Q) J 


3p 


L '^p-3a (n) 


IV^I 


La(n) 


<c(l + ||e|| 


3G^+“\ 

LP{n) J 


p+3<y. 

I 3(p-l) 

\lp{Q.) 


( 21 ) 


where the condition < 7+0 is less restrictive since it requires only 5 + 3 q; < 3(7+0) ^ 

7 > |. The other terms can be estimated similarly so we get taking maximal possible value of 
a = 2(7 — 3 ) that 


II^ 1 |Il 37 - 6 (Q) < C. 


Using this in the already derived estimates for v, c and ^ yields the result of Lemma [TJ 


□ 


Now, we will show that for 7 > 6 we can expect principally better regularity of the solutions, 
this is connected to the fact that in this case we can take according to (|20l) a > 7, so p{q, c) G L^(n), 
for some s > 2. 


Lemma 3 . For 7 > 6 we have for solutions to ( 0 )-([ 3 ) for any 1 < p < +00 


ll’^llivi.p(n) + ll^’llL°°(n) + ll^c||^oo(Q) + ||+|lz,oo(f^) + ||£iL'(c)||^<^^^^ < Cp, 
and c G [0, 1 ] a.e. in U. 

Proof. First, since 7 > 6 we observe certain smoothing effect of (|4]) and ([5]). In what follows, we will 
repeatedly use Holder’s inequality in the third equation and Lemma[2j Indeed, since gv ■ Vc = —p, 

„ _fi fi 67-12 37-6 

and Q G °, v G L° and Vc ^ L we get p & L , ( g^_g + ^ + g ^_^2 > 2 )’ applying 

67-12 

Lemma[2]Vc G which can be again plugged into the third equation in order to get p ^ L r+e ^ 

67-12 

and again Vc G L 1^-7 ^ at least for 6 < 7 < 12, etc. This procedure can be repeated until Vc G L°°, 
since there exists no reasonable solution to the following system of algebraic equations, where P 
corresponds to the expected integrability of Vc, Q to /r and corresponds to gp 


1 11 1 J7 3g(7 - 2 ) 

37 — 6 6 P Q’ 37 — 6 — Qj — Q 


( 22 ) 
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67-12 


So we get that P —>■ +oo, and Q —)• 

67-12 

steps we have Vc E L°°(0), fj, G L -y 


> 3, id est for any 7 > 6 after finite number of such 

67-12 

(O), and Ac E L ^+2 Thus, to summarize 


ll’^lliyi'2(n) + ll£’llL37-6(n) + l|Vc||^oo(Q) + 11 ^ 11 ^ 67 ^^^^ + ||Ac||^ 62 ^^^^ + ||pT'(c)||^ 62 ^^^^ < C. 

^From the last norm we can deduce c E [0,1] a.e. on {g > 0}. On the other hand, c is continuous 
in n and on the set {p = 0} it satisfies the Laplace equation, and therefore maximum principle. 
ThusH c E [0,1] a.e. in O. 

Now, we will use the fact that we work with slip boundary condition, and thus we can deduce 
from the momentum equation the following relation for vorticity (w = curlv), see m 


— uAu = — curl(pv • Vv) — curl(AcVc) + curl(^F) in n, 

k 

oj ■ Ti = -{ 2 x 2 -)v • T 2 on dQ, 

V 

k 

U) ■ T 2 = (2xi -)v • T 2 on dO., 

V 

div w = 0 on dkl, 


where Xn are the curvatures related to the vectors r^,. Note, that n • (Vc® Vc) • = 0 on dkl, 

since Vc • = 0. Now, we would like to show that || < C, for some p > 1. The analysis 

of the above system (the operator defined on the left hand side) can be found in [25]. 

First, we have v • E W 2 ’'^{dQ), so we control a; on the boundary. Further, pF E L^'^~^{Q), 
pv • Vv E LP{Q,), for any p E (l, finally 

II AcVc|| 67-12 < ||Ac|| 67-12 ||Vc||j;^oo(f 2 ) < C. 

L T'+-^ (U) L 7+2 (O) ^ ' 

According to the fact that 1 < < ^22^ even for all 7 > |, we get 

w E IF^’^(n). 

Now, we will proceed in the same manner as in [22|, using the Helmholtz decomposition of the 
velocity field v = curl A + Vcp, define G = — ( 2 z 2 -|- p)A(j) + p{g, c), observe 

VG = —pv • Vv -|- z/Acurl A + (pF -|- AcVc) 
and show that G E kF^’ 27-3 (H) L r-i (H). This further yields 

p7 E L^2/5 (j7), ^ ^ L^^{n) 

and after more iterations p E L°°(H), v E ly^’^(H), c E ly^’^(H), for any 1 < p < 00 as well. □ 

®More precisely, since c is continuous, the set U = {c ^ [0,1]} is open. Considering any ball B{r,xo) C U we get 
that £1 = 0 a.e. in B{r,xo), hence c satishes Ac = 0 in B{r,xo) and can not reach neither maximum, nor minimum 
within this ball. 














3 Approximation 


In this section we define a problem approximating the original one and prove the existence of the 
corresponding solutions. We introduce /i = M/|n|, e>0, a smooth cut-off function K{q) 


K{q) 


1, for Q < m — 1 

< G (0,1), for Q G {m — l,m) 
_0, for g > m, 


a ’’regularized logarithm” which is a function 4 G C^([0, oo)) which is bounded from below by 
log(-v/e) — 1 (t > 1 will be specified later) and 



k{s) 


log(s), 

convex, non-decreasing 


for s > -y/e, 
for s < -y/e, 


0 < -|- s/g (s)) < C, 

0 <s(2/'(s) -|-s/"(s)) < C for a.a. s G [0, oo), 


(23) 


where C is independent of e; further we denote the approximated free energy 


fe{Q, c) = r(£>) -h (ac -h d)le{Q) + Le{c) + b{c), 

C 

where we define Le(c) = f le{s) — 4(1 — s)ds for c G [0,1], and then extend it to whole M as a convex 

0 

function with < — Clog-^, r(^) = and approximated pressure 

Peig, c) = Pb{g) + {ac + d) f K{s)ds{s‘^l{{s)) ds, 

Jo 

where Pb{g) = Jq 'yP~^K(s) ds. 

Our approximation problem then reads 


eg + div(K(g)gv) = eAg -|- eK{g)h, 


(24) 


- div{K{g)gv (8) v) -|- -K{g)gv ■ Vv — zzAv — {ly + r/)V(divv) -|- Vpe{g) 


= K{g)gF + div(Vc ® Vc - - aVc sK'{s)ds, (25) 

2 ./n 


K{g)gv • Vc = - 


dfe 


K{g)gpL = -Ac+K{g)g'^J^ + egK{g)L'{c). 

oc 

Moreover, we supply the equations with additional boundary condition 

• n = 0. 


(26) 

(27) 

(28) 


"We can get such an function e.g. by replacing the logarithm by a suitable affine function for small arguments. 
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Proposition 1. Let e > 0. Suppose that the assumptions of TheoremU\ are satisfied, then there 
exists at least one solution v ^,q to the system ® - ^ 2 ^ with Moreover, we 

have with 1 < q < +oo the following estimates independent of e 


^ II^£'IIl 2 (q) + £ {K{q)q)^^^L'( c) < C{m), 


\\P^\\L'2{n) + II/^IIl2(q) + ||Vc||^^^^^ + ||ve||^yi,2(Q) + 


7 + 1 


{K{ 0 )q) L'fic) 


2j 

Lt+i(O) 


< C. 


(29) 


Proof. The existence of solutions for the approximative system will be deduced by means of the 
following variant of the Leray-Schauder fixed point theorem, see e.g. [ 8 ]. 

Schauder fixed point theorem. Let X be a Banach space, and T a continuous and compact 
mapping T : X X, such that the possible fixed points x = tTx, 0 < t < 1 are bounded in X. 
Then T possesses a fixed point. 

We define spaces 


Mg={w€ • n = 0 on , 

Mg ={me 1T^’'^(0), Vm • n = 0 on dQ} , 

Ng={z£ Vz • n = 0 on , 

and search for g € Mp, v G Mp, and c G Np, with 1 < p < oo, see 0127] for similar considerations 
for the Navier-Stokes system. Let us first concentrate on the continuity equation. 

Lemma 4. The solution operator 5i(v) = g of the problem 


eg + div(iL(p)pv) = eXg + eK(g)h in Q, 

Vp • n = 0 on dQ 

is forp > 3 a continuous compact operator from Mp to W'^’P. Moreover, g > 0, f g dx < M, and 

o 

WqWvV^.p < C'(m,e)(l + l + ||i4/i,p(Q))- 

Proof. The proof can be found in [23|, see also [191122] . We recall here only the idea how to obtain 
the estimates. First, considering the subset {p < 0} C 12 we get £» > 0 a.e. in 12, then integrating 
the approximate continuity equation over 12 yields 

J gdx ~ ^ J 

n n 


so K{g) = 1 a.e. in 12. Further, testing the continuity equation by g yields 



^(s)sds)da; < C\\K{g)g^\\p^^^^.^ l|Vv||^2(t5) , 
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since the last term on the left hand side can be easily bounded we get 


(30) 

□ 




Similarly, for the last two equations we have 
Lemma 5. The solution operator 52 (v) = c of the problem 

-^i = K{q)qS/c-w, 

d f 

—/S.C +K{Q)QeL'(c) = K{q)qpl — K{q)q-;:^ in where g = Si(v), 

oc 

Vc • n = 0 on dTl 

is forp > 3 a well-defined compact operator from Mp to Np. Moreover, c E [0,1], and 

\\^(s)Q\\Lo{n) + !)• 

Proof. The proof is quite similar to the previous Lemma, for the estimates we proceed analogously 
as in the proof of Lemma [2l For constructing the solution we use again Schauder fixed point 
theorem. We consider for fixed g E Mp the mapping c ^ z defined as a solution operator to the 
problem 


—/U = K{g)gVc ■ v, 

-Az + K{g)geL'{z) = K{g)gfi - K{g)g^{g, c), 
Vz • n = 0 on dQ. 


(31) 


The second equation is for e > 0 strictly elliptic, furthermore, its right hand side belongs to 
VF^’^’(n) ■H' L°°{Q), in particular we deduce that K{g)gL'{c) E L°°{Vl) and Vc E VF^’°°(n), hence 
since g and 2 ; are continuous we find K{g)gL"{z) E L°°{Q), put the corresponding term to the right 
hand side, observe that in fact K{g)gfi — K{g)g^ — K{g)geL'{z) E VF^’P(n), and so 2 ; E W^’P{Q), 
and the corresponding mapping is compact by the same reasons as above. In order to get the 
desired estimate independent of e, we test the equation —Ac + K{g)g{^L'fic) + eL'(c)) = K{g)gpi — 
K{g)g(h'{c) + alfig)) by F{L'fic) + eV{c)), where F is increasing function such that CF{C) ~ 

\F\ ICI^ , for some /3 > 0. Note that L'fic) + tL'{c) is non-decreasing function, so we get 


{K{g)gfFP+^)L'{c) 


/ 3+1 
L/3+i(0) 


+ 


{K{g)gfFd+^)L'fic) 


/ 3+1 

L/5+i(r2) 


< 


< 


j\Vc\^F'{L'fic) + eL'{c)){L'f{c) + eL''{c))+K{g)g{L'fic) + eL'{c))F{L'fic) + eL'{c))dx 


K{g)gpF(yL'fic) + eL'{c))dx 


n 


+ 


J K{g)g[b'{c) + alfig))F^L'fic) + eL'{c))dx 


(32) 
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hence 


{K{q)q) p+^L'{c) 


LP+^(Q.) 


+ 


{K{q)q) '3+iL'(c) 


LP+^{n) 

<C{M^,^^^\\Kig)g\\ +1) (33) 


'LP-(/3+1) (Q) 

and using the classical elliptic estimates we obtain as in Lemma [2] 

11 ^ 11 ^ 2 ,^ \\K{g)g\\L,^n) + 1 ), 

especially for p = 2, 


l|c||^2,^ + l|Vc||^^^^^ < C(||p||^2(f^) \\K{g)g\\^,^^^ + 1). 

The issue of existence of the solutions to (j31ll requires some comments. The function L'{-) is singular 
and it keeps the value of z in the interval [0,1]. Thus, we approximate (j31|] by its regularization 
substituting L'(-) by L^(-) which is obtained in the same manner as for in (I23p . The estimates 
are the same, there is no problem to pass to the limit <5 —>■ 0, since we control the second derivatives 
of z. Hence, we ensure that z G [0,1] as well. □ 


Finally, we define the solution operator T : Mp —>■ Mp, T(v) = w of the problem 

— uAw — {u + 7/)V(divw) = — i div(iF(^))^v ® v) — ^K{g)gv ■ Vv — VPe{g) + K{g)gF 

iVcP 

+ div(Vc(8)Vc-— aVc / siF'(s)ds, (34) 

2 Jo 

where g = 5i(v), and c = 52 (v) 


equipped with the boundary condition 

w • n = 0, n • T(w) • Tn + k'w • r = 0, on dO.. 

Lemma 6. P is continuous and compact operator from Mp to Mp for p > 3. 

Idea of the proof. It is again one more time strictly elliptic system with right hand side which 
belongs at least to the L'P{Vt) and it contains at most first order derivatives of g, v and at most 
second order derivatives of c, see e.g. m for similar considerations. □ 

Finally, we will verify that all possible solutions of tT(v) = v, for t G [0,1] are bounded in Mp 
independently of t. Testing the momentum equation by v yields (with usage of the last equation 
tested by V • Vc) 


t J v-Vpe(^, c) dx + J jVvj^ + (i 2 + T/) jdivvj^dx + J A:(|vi • ti|^ + |v 2 • T2|^)d5 

n n dfi 

— K{g)gF-vdx + t J ^K(g)gfj,Vc ■ V — K(g)g^^Vc ■ v^dx 


n 


n 


d- ^ J (^—egK{g)L'{c)Vc-v — aVc-v J siF'(s)ds^dx, (35) 


n 
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from the equation for concentration we get f /u^dx = — f K{g)gfiv ■ Vcdx. Next, 

n n 


J v ■'Vpe{g,c) dx = jw-Vgi^g'^ ^ K{g) + {ac + d)K{g)-^[g^l'^{g))^dx 
n n 

+ J aVc-v j iir(s)^(s^/'(s))ds dx (36) 

Q. ° 


and 

j K{g)gv ■Vc^{g,c)dx = J K{g)gv ■ V {fe + g^)dx - J K{g)g^al[{g)w ■ Vc dx 

ft ft n 

— J K{g)gv ■ Vg(2T'{g) + gT''{g) + (ac + d)l'^{g) + (ac + d)-^{gl'^{g))'^dx 
n 

= J div(^K{g)gv{f^ + g^)^dx - J dW{K{g)gv) (/^ + £'^)dx ( 37 ) 

ft ft 

— J aVc • V J —(it'(s)s^Z((s))ds dx 

ft ° 

/ I 

K{g)'v ■ Vg—{g'^r'(g) + g^l[{g){ac + fi))dx. 

ft 

The first term on the right hand side of (1371) can be eliminated by boundary conditions, for the 
second one we will use the continuity equation (tested by (/e + Q^))- The third one cancels out 
with the last terms of (1551) and (f55D , and the last term is exactly the same as the main part of (1551) . 
Thus, summing up the resulting inequalities with appropriate powers of t we get 


h\\wf^{ft)+t\\^^\\\■^(ft)+'t J (^Q{fe + g^)dx <t J K{g)gF -vdx -t J egK{g)L'{c)^r -Vcdx 

+ t J {eAg +eK{g)h){fe +g^)dx, (38) 

ft 

where we have introduced the notation Qk = {K{q)q > 0}nO. We estimate the e-dependent terms. 
K{Q)QL'ic)'Vc ■ vdx = — J eK {g) gv ■ V [L{c)) dx 

ftK flK 

J €gK{g)v ■ nL{c)dS + J ediv{K{g)gv)L{c)dx 


(39) 


dVtK 

<V^\\^Q\\L^ft) \\QK'{g)+K{g)\\^,^^^^ I|v||i6(n) ^ll^(c)||i.^(nj^) 


+ \\K{g)g\\L2(^Q-^ l|Vv||^2(Q) ^ ||-t'(c)||ioo(Q^) 


, 3/4 


and we can use ||T(c)||£cx)(Qj^) < C{-^\\{K{g)g)^/^f ^~^^^+ 1) with estimate ([M])- 
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Concerning the term eK{g){f^ + Q^) we have three parts. First, re(^)) + ^3r'(^3), which can be 
bounded by the corresponding term on the left hand side, see e{gfe)- Second, {ac + d)(l + le{g)), 
which has good sign on the set {£» < e“^} and is bounded hy c{l + g)K{g) on {^ > e“^}. The third 
term e(^L^{c) + b{c)) can be bounded according to our definition of For Ag{f^ + ^?^) we have 

J eAg{f^ + g^)dx =-e j\Vg\^ {yg^~‘^ + {ac +d){gl”{g)+ 2l[{g))^dx 
n n 

-e J Vg-Vc{l'^{g)g + b'{c) + L'^ic))dx. (40) 

Q 

The first integral has a good sign. Indeed, g'^~'^ as well as gl'l{g) + 2 ^((^ 3 ) are non-negative (for large 
arguments we have gl”{g) + 2l'^{g) = -i > 0 , in the other case the conclusion is obtained from the fact 
that for small arguments is increasing and convex, see (l23l) i. Note also that we have (l30]l and for 
e > 0 the approximated version of (l5|) yields c G [0,1] as soon as we control 11.IF(^))^)||^ 27 (^ 2 ) ll/^llL 2 (f 2 ). 
For the rest we get denoting 14(c) = b'{c) + L'^{c) and U = ||v||^i, 2 (f^) +111/^1112(^2) that 


t j e|V£-| |Vc| \l'^{g)g + Veic)\dx < q\\l 2 (^^) l|Vc||i 2 (f 2 ) v^(||le(e')£> + K(c)||^^(f^)) 

Q 

< tree \\K{g)g^^/,]^^ mQ)Q\\L.in) Mlhu) 

< tre ||iF(^^)^lll(^) \\K{g)g\r^^) (41) 

where we have used the choice of 4 which guarantees that ^ ||I((^)i?|lx,oo(Q) + || 14 (c)|| 2 ^oo(q)) < C 

independently of e. Thus, 


[/ < C (1 + \\Kig)g\\l,;,^^^ + e\\K{g)g\r^^r'^^ + \\K{g)g\\L3^^r'^^)' 

However ||iF(^)^)||^6(Q) is definitely finite so finally, 

e < -I- ||^F(£))^>||l6/5(Q)), < C{\\K{g)g\\j^g(^Q^ \\f'\\L'2{n) + !)■ 

The Bogowskii estimates go along exactly the same lines as in the a priori approach as soon as we 
observe that 


/ 


av • Vc f tK'{t)dt^dx 

Jo 






(42) 


so 


I i?e 112,37-6(0) < C, 


independently of m and e. Furthermore, by the same iteration process applied on the last two 
equations of ©-([SI) as in the a priori approach we can deduce that 




ell 67 + 

L^{n) 


IVC, 




T 


67 ^ 

L^{n) 


C. 


(43) 
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Having these estimates in hands and noting that \\K{g)g\\j^ao^Q'^ < C{m) we can apply the elliptic 
theory on the equation (j34|] and get the estimate of fixed points of T in Mp. This completes the 
proof of Proposition [H □ 


4 Artificial diffusion limit 


The last section of the paper is dedicated to the proof of convergence of the constructed approxi¬ 
mative solutions to a weak solution to the original system. As usual the key part is related to the 
proof of the strong convergence of the densities. 

Thanks to the estimates (f2^ . ()43]) we extract from the family v^, /i^, Cg) subsequences which 
converge in the corresponding spaces as e —0 + . Namelyjl 




Qe 

K{ge)Qe 

tK'{t)dt 


eK{g^)g^L'{ce) 


V in 
g in 

la^ in 


tK'{t)dt in L°°{Q), 


Jo 

cin 

eK{g)gL'{c) in 
g in 


V, ^ V in 

Pei0,c) -^*p{g,c) in L°°{n), 
K{g,) 

hM ^*^in 

Vc, ^ Vc in 

K{g,)g,L[{c,) - A®L'(c) in 
K{ge)gePe ^ K{g)gfi in 


Thus, we get 


div(iir(^))^v) =0, 


- div(A'(^>)gv ® v) -|- -K{g)g v • Vv — i/Av — [v + r/)V(divv) + Vpe{g, c) 


_ iVcl 

= K{g)gY + div(Vc (g) Vc-/ tK'{t)dt, 

2 Jo 


K{g)gv- Vc = - p, 


K{g)gp = - Ac + h^{g) + K{g)g[L'{c) + b'{c)) + eK{g)gL'{c), 
V • n = 0, Vc • n = 0, 
n • T(c, v) ■ Tn + kv ■ Tn = 0, on dQ, 


where /ie(s) = as ■ le{s). Note that, due to the high regularity we have the pointwise convergence 
of concentrations. In order to show the pointwise convergence of densities we need to investigate 
the momentum equation, especially its potential part defining the effective viscous flux. Let us 
decompose the velocity field v using Helmholtz decomposition, i.e. 


V = V<y9 -|- curl A, 

®We denote a weak limit of nonlinear expressions {we(pe, Ve, pe, Ce)} by «;(£>, v, p, c). 
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where 


curl curl A = curl v = oj in 

div curl A = 0 in n, and 

curl A • n = 0 on 917, 


A(/j = curlv in 17, 
V(^ • n = 0 on 917, 


J <y9 dx = 0. 

n 


For the stream fnnction A we have good estimates 

||VcurlA||^,(f^) < IIV^cnrlAll^,^^^ < , 


and since to solves 


— uAu) = — cnrl(A(^)^>v • Vv) + curl(A(^))^3F) + curl^AcVc — aVc J tK'{t)dtJ , (44) 

we have ||<^||vi/i. 9 (q) < C{m). Similarly we also decompose the approximative velocity field as 
Vg = + curl Ag and deduce due to the slip boundary conditions the following problem for 

vorticity 


-uAcOe =cuiI(^K{q^)q^F - ■ Vv^ - 

1 

+e-Qe^e + Ac^VCe — aVCe / tK' 
2 Jo 


=:Hi 

k 

w, • n = -(2X2-)ve • T 2 on 917, 

V 

k 

uj^-T 2 = (2x1 - -)v, • T 2 on 917, 
divwe = 0 on 917. 


- curl(-eA£»,Ve) =: Hi + H2, 


= :H2 


The structure of the right hand side of the first equation enables us to consider = a;)* + + w, 

as a sum of solutions to three particular systems, namely 


uAu}^ =0, 

• Ti = - (2x2 - -)ve 
V 

ujI-T2=- (2X1 - ^)Ve 
diva;° =0, 


vAu)\ =Hi, 

T2, OjI ■ Tl =0, 

Tl, 0 j\ • T2 =0, 

divw^ =0, 


vAuJl = H 2 in 17, 
cjg • Tl = 0 on 917, 

cjg • T 2 = 0 on 917, 
divcjg = 0 on 917. (45) 


Lemma 7. If the vorticity = 0 ;° + cjg + uj‘1 solves then we have 






\Li{n) 


+ 1 ^), for q G 

< C'(m)e^'^^, for q£ [1,2]. 


2 , 


67 

6 + 7 


(46) 

(47) 
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Proof. Following closely the corresponding considerations in [23], we deduce that (see 
li^c ||vKl’9(0) — ^ ll^e < C* ||Ve||y[/i,9(Q) , 

but according to (l43]l for any q such that 2 < q < we have 

||Ve||wl,.(Q) < Ci\\PeiQe,Ce)\\L^Q) + || Ac,Vc, + 1) < C {l + , 

SO we concentrate on col. If we denote r = > 3 we get by interpolation 

|L,1|| ^ ("] I \\ T/' ( ^ \ ^ V7-.T II 

ll^e |liyi,g(Q) 


<C(1 

+ 

\\K{ge)geVe'^Ve\\Li{n) + l|Ac,Vce 

L9(0)) 


<C(1 

+ 

m 

l|v6|lLcx,(f^) ||Vv, 

E Ll9(0)) 







2r-6 

r 

2(r- 

-9) 

'■(9-2) 

<C(l 

+ 

m 

i|v.iii'.(^, iivv. 

II II Vv II 

IIl'-(O) II '^^dl 

9(e- 2) ||\^„ 
I,2(n) II VV 

II 9('--2) 

eIIl’'(o) 

<C(l 

+ 

m 

r ,r(q- 

II Vv ||3>— 

< C(1 + 

] 

m 

1 1 4rQ —6r 
'-+39(r-2) 

7(1-^)) 


(48) 


Finally, for the last part we get for q < 2 

Now, we are approaching the key definition of the effective viscous flux. Inserting the Helmholtz 
decomposition into the approximative momentum equation yields 


V{-{2u + ri)A(p^+p^{g^,Ce)) = z/AcurlA + K{g^)Q^F - K{g^)g^w^ ■ Vv, - ^eA^.v^ 

1 1 iVc P 

- -ehK{ge)ve + e-ge^e + div(Vce ® Vce-- CL^Ce / tK'{t)dt, (49) 

2 2 2 Jq 

and we introduce the fundamental quantity 

G, = -{2u+ r])A(p+p^{g^,Ce) = -{2u + rj) div + p^{g^, Ce). (50) 

Similarly, inserting the Helmholtz decomposition into the limit momentum equation we obtain 
(with usage of the fact that due to the continuity equation we have K{g)gVw = K{g)gv ■ Vv) 
that 


V (—(2z^ + p)Aip + p{g, c) ) = i/Acurl A + A(^)^>F — K{g)gv ■ Vv 

iVcl^ ~T^ 

+ div(Vc(8)Vc- ^^tj—aVc tK'{t)dt, (51) 

2 Jo 

hence we define the limit version of effective viscous flux by 


G = -{2u + r])A(p +p{g,c) = -{2u + rj) divv +p{g, c). (52) 


Note we have control of f G^dx = f pe{ge-,Ce) dx and f Gdx = f p(g,c) dx. Further we state the 

no o o 

most important features of the effective viscous flux. 
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Lemma 8. There exists a subsequence such that 


Ge ^ G (strongly) in 


and 


||G||^oo(f^) < C'(C)(l +m^+ a'+C), for any C S [o, 


7 — 6 


(53) 


Proof. Let us decompose G^ to G^ = G) + G'(, where G^ contains the ’’strongly e-dependent” terms 
of the right hand side of (l50]) . namely VG^ = — e^A^^Ve — z^curlwg with f G^dx = 0, so 

n 

ll‘^e|L 9 (Q) < +z^||curlw2||^^_^^^^P < C(m)e^/^ for q G [1,2], 

Using once more the estimates from Lemma [71 we observe that 1/^ G^dx] < G and 


iGil 


w^’i(n) 


< G(1 -|-for g 


G 


67 

6 - 1-7 


Therefore, since 7 > 6 we have, at least for a suitably chosen subsequence 


(54) 


G 


G^ (strongly) in L°°(Q) and G^ ^ 0 (strongly) in L^(fl). 


Thus, 

Gf; = G) + G^ —i- G^ = G (strongly) in for q G [1, 2], 

Finally, setting (7 = 3-1- 2 ^^ (l54]i we get the desired conclusion 

||G||^^(^) <G((7)||G||^i,g(^) <G(C)(l + mi+^+<). □ 

Now, we are ready to show that we are able to choose m in such a way that actually K{q)q = g 
a.e. in n. This will be an immediate consequence of the following lemma. 

Lemma 9. There exists niQ such that 

iryi _ 

-(m - 3)"^ - ||G||^oo(m >1 form> mg, (55) 

and at least for a subsequence lim > m — 3}| = 0. 

€->■0 

Proof. Let us define a smooth non-increasing function N : [0, 00) —>■ [0,1] such that 


N{t) 


1 for t G [0, m — 3], 

< G [0,1] for t G (m — 3, m — 2), 

0 for t G [m — 2, 00 ), 


and multiply the approximative continuity equation by N\g^), for some suitable power / G N in 
order to get after some manipulations 


Be 


tlN^ ^(t)A'(t)dtdivvgdx > iie 


n 0 


(56) 


18 











with Rf; = ef N^(g^)Agf^ + (h — g^)N^(gf;)dx, -R^ — )• 0 as e ^ 0 @ Further by definition of Gg we 
have 


Q 


Be 

-(m-3) j lN'‘~^{t)N'{t)d-^pe{ge,Ce)d. 


X < m 


and 


n 0 


m — 3 


Be 


-lN^-\t)N'{t)dt\ G,dx 


n 0 


+ Re, (57) 


m 


{l - N'-{g^))p^{g^,Ce)dx < / (l - iV'(£i£)) |G,| dx + |R,| . (58) 


{ee>m-3} {ee>m-3} 

Recalling the structure of the pressure we have according to (f23]l and the fact that € [0,1] 

'•Be 


Pe{Qe,Ce) = PbiOe) + (aC^ + d) K dt > Pb{ge), 


(59) 


'0 


which yields 
m — 3 


m — 3 


m 


{m-3y\{ge >m-3}\ - — \\Pe{Qe,Ce)\\L 2 ^^) N^{ge) 


L2(n) 


^ IIG'IIloo(q) !{£>£ > m - 3}| + ||G - + \Re\ . (60) 

According to (f53]l we are able to choose mo satisfying (f55]l . yielding 


|{£>, > m-3}\ <G{ N^{ge) 


P{{Be>'m—3}) 


+ IIG - G'ellil(n) + l-^el)- 


(61) 


However, the last two terms tend to zero as e —0+ and having fixed e > 0 ||-^*(^^e)||^ 2 ({g 
tends to zero as I —>■ +oo as well. Thus, Lemma [9] is proved. □ 

Finally, we deduce the pointwise convergence of the densities. Our main aim is to show that 


Pbi0)0 = Pb{0)0, 

which will further lead to g^ ^ g strongly in L^(H) for q < oo. 
Lemma 10. The weak limits satisfy 

J p{Q,c)gdx < J Ggdx. 


(62) 


Proof. Testing the approximative continuity equation by logm — log(^£ + 5) and taking the limit 
for 5 —)• 0 +, we get (see [ 22] 1 

/ A'(^e)ve • V^edx > eC'(m) 


Note that —d f N’' ^{q^N' (qc) \V gdCA x > 0 
n 
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and by Lemma [9] f divv^dx < with —)• 0 as e ^ 0 + . Hence the definition of G yields 


Ve{Qe, Ce)£'edx < J G^Q^dx - {2v + r/)R, 


and passing to the limit with e we get (l62]l , since according to Lemma [8] Gg = Gg. 


□ 


Lemma 11. 

J p(g,c)gdx = J Ggdx. (63) 

ft 

Proof. First, with the usage of Fridrichs’ commutator lemma we are able to approximate g by 
smooth bounded functions gn and deduce that 

J £»divv + V • Vfjdx = 0. 

n 

Further testing continuity equation by log(^3ri + <5) — log 5 and then passing to the limit at first with 

n —>■ oo and then with 5 —>• 0+ gives us f v ■ V^^dx = 0, hence f ^divvdx = 0 as well. Thus, using 

n n 

again the properties of G we can conclude (l63]l . 

□ 


Further, the strong convergence of c and the convexity of s —)• and s —)• s‘^l[{s) gives us 

p{q, c)g < p{g,c)g, 

which combined with (j62l) and (16311 yields 

p{g,c)g = p{g,c)g a.e. in 
Thus, ^'^+1 = gg'^ and g^ ^ g strongly in 

Now, we move our attention to the last two equations of system (I2|)-(l5|) and show that due to 
strong convergence g^ and all the remaining nonlinearities can be identified, so we have indeed 
obtained the solutions to our original system. This completes the proof of our main theorem. 

Remark 2. In order to prove the existence of only weak solutions of the system under consideration 
for 6 > 7 > 3 we can use the following idea. We modify the pressure psig, c) = p{g, c) + with 
F > 6, so we can apply the already proven result and then using the a priori estimates derived in 
the heuristic approach pass to the limit with <5^0 + . This limit passage can be performed in 
the same spirit as in the case of the Navier-Stokes system, using nowadays ’’standard” techniques 
due to Lions, and Feireisl, see e.g. [26] with 9 replaced by c. The compactness of the additional 
stress in the momentum equation and in the additional equations is just easy application of the 
Rellich-Kondrachov compactness theorem due to the uniform bound of Ac. 
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